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We analyze spatio-temporal structures in the Gross-Pitaevskii equation to study the dynamics of 
quasi-one-dimensional Bose-Einstein condensates (BECs) with mean-field interactions. A coherent 
structure ansatz yields a parametrically forced nonlinear oscillator, to which we apply Lindstedt's 
method and multiple-scale perturbation theory to determine the dependence of the intensity of 
periodic orbits ("modulated amplitude waves") on their wave number. We explore BEC band 
structure in detail using Hamiltonian perturbation theory and supporting numerical simulations. 
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At low temperatures, particles in a gas can reside in the 
same quantu m (g round ) state, forming a Bose-Einstein 
condensate Q, OH liii UJl UM, Lili • When considering only 
two-body, mean- field interactions, the condensate wave- 
function -0 satisfies the Gross-Pitaevskii (GP) equation, 
a cubic nonlinear Schrodinger equation (NLS) with an 
external potential: 

ihip t = -[h 2 /(2m)]i> xx + g\iP\ 2 ip + V(x)il> , (1) 

where m is the mass of a gas particle, V(x) is an external 
potential, g = [47r/i 2 a/m][l +0(l 2 )], a is the (two-body) 
s-wave scattering length and i = -\/|?/>| 2 |a| 3 is the dilute 
gas parameter I23. The quantity a is determined 
by the atomic species in the condensate. Interactions 
between atoms are repulsive when a > and attractive 
when a < 0. When a » 0, one is in the ideal gas regime. 

The quasi-one-dimensional (quasi-ld) regime em- 
ployed in JQ) is suitable when the transverse dimensions 
of the condensate are on the order of its healing length 
and its longitudinal dimension is much larger than its 
transverse ones |6|, tZl ISl LLiJ ■ ^ n this situation, one em- 
ploys the Id limit of a 3d mean-field theory rather than a 
true Id mean-field theory, as would be appropriate were 
the tranverse dimension on the order of the atomic inter- 
action length or the atomic size. 

In this paper, we examine uniformly propagating co- 
herent structures by applying the ansatz iji(x — vt, t) — 
R(x — vt) exp [i [9(x — vt) — cut}), where R = |^>| is the 
magnitude (intensity) of the wavefunction, v is the veloc- 
ity of the coherent structure, 9{x) determines its phase, 
and lo is the temporal frequency (chemical potential). 
Considering a coordinate system that travels with speed 
v (by defining x' = x — vt and relabeling X clS X ) yields 



t/j(x, t) — R(x) exp (i [9{x) — cut}) . 



(2) 



When the (temporally periodic) coherent structure is 
also spatially periodic, it is called a modulated amplitude 
wave (MAW) . The orbital stability of coherent structures 
(0) for the GP with elliptic potentials has been studied by 
Bronski and co-authors [a, 13 ■ To obtain information 
about sinusoidal potentials, one takes the limit as the el- 
liptic modulus k approaches zero [2^,|33- When V(x) is 
periodic, the resulting MAWs generalize the Bloch modes 
that occur in linear systems with periodic potentials, 
as one is considering a nonlinear Floquet-Bloch theory 
rather than a linear one jj, Oj [T3, 12a, yj] . 

In this paper, we employ phase space methods and 
Hamiltonian perturbation theory to examine the band 
structure of such MAWs. Prior work in this area has 
utilized numerical simulations 0, 0, ■ 

The novelty of our work lies in its illumination of BEC 
band structure through the use of perturbation theory 
and supporting numerical simulations to examine 2n : 1 
spatial subharmonic resonances in BECs in period lat- 
tices. Such resonances correspond to spatially periodic 
solutions of period 2n and generalize the 'period doubled' 
states studied by Machholm, et. a /. 12911 and observed ex- 
perimentally by Cataliotti, et. qL.|l2j| 

Inserting J2J into the GP QJ, equating real and imag- 
inary parts, and simplifying yields 



R — S , 

c 2 2mujR 2mq 

<j' — 1 

i? 3 h h 2 



R" 



(3) 
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The parameter c is defined via the relation 0'{x) = c/i? 2 , 
which is an expression of conservation of angular mo- 
mentum Null angular momentum solutions, which 
constitute an important special case, satisfy c = 0. 

When V(x) = 0, the two-dimensional dynamical sys- 
tem © is autonomous and hence integrable. Its equi- 
libria and the stability thereof are discussed in detail 
in Ref. 1^3]. When g > 0, lj > 0, and c = 0, which 
is the primary case for which we study the band struc- 
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situations, such as for 85 Rb (for which a — —0.9), are 
examined in Ref. [3^ . 




FIG. 1: Phase portrait of a repulsive BEC with no external 
potential and w = 10. In this plot, the two-body scatter- 
ing length is a = 0.072 nm, obtained for atomic hydrogen 
C~H) |20|. Orbits inside the separatrix have bounded intensity 
R(x) with increasing period as one approaches the separatrix. 
The quantities R and S — R' are scaled quantities obtained 
with m = 1/2 and h = 1. 



ture, one obtains a neutrally stable equilibrium (a cen- 
ter) at (R, S) — (0, 0) and unstable equilibria (saddles) 
at (±y/fko/g,0). See Fig. □ 

We employed Lindstedt's method to study the depen- 
dence of the wave number of periodic orbits (centered at 
the origin) of 10 on the intensity R when V(x) =0 |34j . 
We assumed g = eg, where e < 1 and g = 0(1). The 
wave number is then a = 1 — 3gA 2 / (8tuh) +0(e 2 ), where 
R(£) = Ro(£) +Q(e), £ := ax, i? (£) = Acos(/3£), and 
(3 := yj2muj/h. 

To study the wave number-intensity relations of pe- 
riodic orbits in the presence of external potentials, we 
expand the spatial variable x in multiple scales. We 
define "stretched space" ^ := ai as in the integrable 
situation and "slow space" r\ :— ex. We consider po- 
tentials of the form V(x) = eV(£,,rj), where V(£,r)) = 
Vo sin [«(£ — £o)] + Vi(rj) and Vi(t]), which is of order 
0(1), is arbitrary but slowly varying. Cases of partic- 
ular interest include Vi(ry) = (periodic potential) and 
Viiv) = Vfiiv ~ Vo) 2 (superposition of periodic and har- 
monic potentials). When Vh -C Vo, this latter potential 
is dominated by its periodic contribution for many peri- 
ods [Til Hsj. The wavenumber parameter is k — 2tt/T, 
where T represents the periodicity of the underlying 
lattice. Optical lattices with more than twenty peri- 
ods have now been created experimentally. Spatially 
periodic potentials, which is the primary case we con- 
sider, have been employed in experimental studies of 
BECs [2y 23 | . They have also been studied theoretically 
in RefsTHE S Hi H E E3 • An example of a coher- 
ent structure for hydrogen in the presence of a periodic 
lattice is depicted in Fig. [3 Coherent structures in other 




FIG. 2: An example of a spatially quasiperiodic coherent 
structure for H in a sinusoidal lattice. As in Fig. |u = 10, 
a = 0.072, c = 0, m = 1/2, and h = 1. Additionally, 
Vo = 10. This plot, which was obtained horn the coher- 
ent structure ansatz, depicts Re(?/>), with the initial point 
(-R(O), 5(0)) = (0.05, 0.05) inside the separatrix in Fig. [I] The 
darkest portions are the most negative, and the lightest are 
the most positive. The quantities R and S are defined as in 
Fig.Q 

When k 7^ ±2/3, the wave number-intensity relation 
for periodic orbits of © is 



a(C) 



1 



2lu 



V^+Oie 2 



(4) 



where Rq(€,v) = A (v) cos(/3£) + B(T})sm(J3g) and C 2 = 
A 2 + B 2 is a constant. (Note that in equation Q and 
in our forthcoming discussion, the small paramater e has 
been absorbed back into the constants, so that we need 
not utilize bars over these quantities.) When k = ±2/3, 
one obtains an extra term due to 2 : 1 resonance: 



a R (C)=a(C)T^ L + 0(e 2 ) 



(5) 



where C is defined as before but is no longer constant, 
and the sign of Vq/Auj alternates depending on which 
equilibrium of the slow dynamics of © one is consid- 
ering [32j. Note that the parameter e has been absorbed 
back into the external potential in equations Q and J5| • 
To examine the band structure of BECs in periodic lat- 
tices, we expand (0 with V{x) = 0, c = in terms of its 
exact elliptic function solutions, convert to action-angle 
variables, and apply several canonical transformations 
to obtain a "resonance Hamiltonian," which we study 
both analytically and numerically. Using elliptic func- 
tions rather than trigonometric functions allows one to 
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analyze 2n : 1 subharmonic resonances with a leading- 
order perturbation expansion l32l l3Sl | . We focus here on 
the case g > 0, u > 0. In Ref. [22], we also discuss the im- 
plication of the work of Zounes and Rand 38] for the case 
g < 0, lo > and briefly consider the technically more 
complicated case g < 0, u < 0. Refs. lEElES] concen- 
trated on numerical studies of band structure. In con- 
trast, we employ Hamiltonian perturbation theory and 
study the band structure of BECs in periodic potentials 
both analytically and numerically. 

Let Co = 7r/(2«) and Vi(x) = 0, so that V{x) = 
(2m/h)Vo cos(kx). Equation is then written R" + 
SR + aR 3 + eRcos(kx) = 0, where S = 2moj/h > 0, 
a = -2mg/h 2 < 0, and e = -(2m/h)V . (Note that the 
perturbation parameter e is not the same as the param- 
eter e employed earlier.) When Vb = 0, one obtains the 
exact elliptic function solution 



R = ppcn(u, k) 



(G) 



where u = uix+u , u\ — S+ap 2 , k 2 — (ap 2 )/(2[S+ap 2 ]), 
ui > 0, p > 0, k 2 € R, and fi G {-1, 1}. The initial 
condition parameter uq can be set to without loss of 
generality. We consider u± S M. in order to study periodic 
solutions inside the separatrix (depicted in Fig. . One 
need not retain the parameter p, to do this, so we set it 
to unity. Because k 2 e [— oo,0], we utilized the recipro- 
cal complementary modulus transformation in deriving 
© UE1 13 H3 Defining X = VSx and r = y^S/^R 
and denoting ' := d/d\, the equations of motion take the 
form 

r" +r-r 3 + eS^ 1 cos (k<5" 1/2 x) r = Q , (7) 
with the corresponding Hamiltonian (using s := r') 



H(r,s,x) = 



1 2 1 4 

-r z r 4 



^r 2 cos ( 4=X ) • (8) 



The frequency of a given periodic orbit is Cl(k) = 
7T\/l — p 2 1 (2K(k)), where K(k) is the complete elliptic 
integral of the first kind [3(|. One thereby obtains the 
action Him |H|33 



J = 



3tt 



E(fc) - ( 1 - y ] K(k) 



(9) 



where E'(fc) is the complete elliptic integral of the second 
kind, and the conjugate angle $ := $(0) + Cl(k)x- The 
frequency O(fc) decreases monotonically as fc 2 goes from 
— oo to (as one goes from the separatrix to (0, 0)). 

After applying several near-identity canonical transfor- 
mations and expanding elliptic functions in Fourier se- 
ries |32| , one obtains an autonomous resonance Hamilto- 
nian K n (Y,£) (in action-angle coordinates) for the 2n:l 
resonance band, 



K n = Y-Y 2 - Kj ^ 



2nVS 



— YB n (Y) cos 
2d 



2n£ 



(10) 



where B n are obtained from Fourier coefficients 
The band associated with 2n : 1 subharmonic spatial res- 
onances is present when k/VS < 2n. 

To obtain an analytical description of these resonance 
bands, we note that such bands emerge from the action 
Y = Y n , the location of the nth resonance torus in phase 
space, which is determined by k/\5 = 2nQ(Y n ). The 
saddles and centers of this resonance band are given, re- 
spectively, by Y s =Y n -\AY\ and Y c =Y n +\AY\, 



AY = *26 



B n {Y n )+Y n B' n {Y n ) 



fi(y„)Vl " 2Y n K'(Y n ) 



1 



(11) 



where K'(Y) := 2K' {k{Y)) / tt , K'(k) := K(y/1 - fc 2 ), 
AY" > when n is even, and AY < when n is odd. 
The width of resonance bands is 



W = 2 



Y n B n (Y n ) 



2nVS 



J"{Y n ) 



1/2 



(12) 



0.32- 




-3-2-10 1 2 3 

FIG. 3: Resonance Hamiltonian K2 for k = 2.5, 5 = 1, and 
e = 0.01. The vertical axis is in units of action Y, and the 
horizontal axis is in units of £/J'(Y). As discussed in the 
text, these quantities are scaled and physically unitless. 

We compare our analytical results with numerical sim- 
ulations in (R,S) coordinates with a = —1, m = 1/2, 
and H = 1. For example, when k — 2.5 and 5=1, the 4:1 
resonances are the lowest-order resonances present. The 
resonance Hamiltonian K2 is depicted in Fig. [3| for e = 
0.01, and the corresponding Poincare section is shown in 
Fig. 0] From (JTDJ, one predicts that the i?-axis saddles 
are located at (R,S) — (±0.86364,0), which is rather 
close to the true value of about (±0.88,0). The S'-axis 
saddles are predicted to occur at (R, S) = (0, ±0.68389), 
whereas the true value is about (0. ±0.687). Numerous 
other examples are studied in Ref. |32| . 
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FIG. 4: Lower right corner of a Poincare section for « = 2.5, 
5=1, and e = 0.01. Note that there is no 2 : 1 resonance band 
for this choice of (k, S). The 4 : 1 resonance is depicted. (Three 
additional copies of this structure appear in the Poincare sec- 
tion.) Recall that the scaled quantities R and S are unitless. 



At the center of the KAM islands, we observe 'period- 
multiplied' states. When n — 1, these correspond to the 
period-doubled states studied recently by Machholm and 



We note, finally, that to analyze three-body interac- 
tions (which is necessary, for example, to examine Fesh- 
bach resonances one has to take dissipative effects 
into account |'25l |35|. In the present paper, we studied 
only two-body interactions. 

In sum, we employed Lindstedt's method and multiple 
scale analysis to establish wave number-intensity rela- 
tions for MAWs of BECs in periodic lattices. With this 
approach, we studied 2 : 1 spatial resonances and illus- 
trated the utility of phase space analysis and perturba- 
tion theory for the study of band structure as well as the 
structure of modulated amplitude waves in BECs. Using 
a more technically demanding perturbative approach re- 
lying on the elliptic function structure of solutions of the 
integrable GP yields the 2n : 1 spatial resonances (band 
structure) of MAWs, which we studied in considerable 
detail both analytically and numerically. 

Valuable conversations with Eric Braaten, Michael 
Chapman, Mark Edwards, Nicolas Gamier, Brian 
Kennedy, Panos Kevrekidis, Yueheng Lan, Boris Mal- 
omed, Igor Mezic, Peter Mucha, and Dan Stamper-Kurn 
are gratefully acknowledged. We are especially grateful 
to Jared Bronski, Richard Rand, and Li You for several 
extensive discourses. 
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